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Abstract. Massive protostars have associated bipolar outflows with velocities 
of hundreds of km s~^. Such outflows produce strong shocks when interact with 
the ambient medium leading to regions of non-thermal radio emission. Under 
certain conditions, the population of relativistic particles accelerated at the ter- 
minal shocks of the protostellar jets can produce significant 7-ray emission. We 
estimate the conditions necessary for high-energy emission in the non-thermal 
hot spots of jets associated with massive protostars embedded in dense molec- 
ular clouds. Our results show that particle-matter interactions can lead to the 
detection of molecular clouds hosting massive young stellar objects by the Fermi 
satellite at MeV-GeV energies and even by Cherenkov telescope arrays in the 
GeV-TeV range. Astronomy at 7-rays can be used to probe the physical condi- 
tions in star forming regions and particle acceleration processes in the complex 
environment of massive molecular clouds. 



1. Introduction 



Massive stars appear to be formed in the dense cores of massive cold clouds 
(Garay &: Lizano 1999, and references therein). The accumulation of gas in the 
core might proceed through previous stages of fragmentation and coalescence 
progressively generating a massive protostar that then accretes from the envi- 
ronment (e.g. Bonnell, Bate & Zinnecker 1997). An alternative picture is direct 
accretion onto a central object of very high mass (e.g. Rodriguez et al. 2008; 
see Shu et al. 1987 for the basic mechanism). In any case, the prestellar core 
is expected to have angular momentum which would lead to the formation of 
an accretion disk. The strong magnetic fields inside the cloud that thread the 
disk should be pulled toward the protostar and twisted by the rotation giving 
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rise to a magnetic tower, with the consequent outflows, as shown by numerical 
simulations (e.g. Banerjee & Pudritz 2007). 

Evidence supporting the existence of molecular outflows is found through 
methanol mascrs, which are likely associated with shocks formed by the interac- 
tion with the external medium (e.g. Plambeck & Menten 1990). However, the 
most important piece of evidence for outflows is in the form of thermal radio jets. 
These jets are observed to propagate through the cloud material along distances 
of a fraction of a parsec in some cases (e.g. Martf et al. 1993). At the end point 
of the jets, hot spots due to the terminal shocks are observed in several sources. 
In a few cases, these hot spots are clearly non-thermal, indicating the presence 
of relativistic electrons that produce synchrotron radiation (e.g. Araudo et al. 
2007, 2008). 

A population of relativistic particles in the complex environment of the 

massive molecular cloud where the protostar is being formed will produce high- 
energy radiation through a variety of processes: inverse Compton (IC) up- 
scattering of IR photons from the cloud, relativistic Bremsstrahlung, and, if 
protons arc accelerated at the shock as well, inelastic proton-proton (pp) colli- 
sions. If such a radiation is detectable, 7-ray astronomy can be used to shed light 
on the star forming process, the protostar environment, and cosmic ray accelera- 
tion inside molecular clouds. In the present contribution we want to discuss the 
conditions under what massive protostars can produce 7 rays detectable by the 
Fermi satellite in the immediate future, and also whether Cherenkov telescopes 
may eventually detect these sources. The models presented here incorporate 
many refinements from the simpler approach published by Araudo et al. (2007) . 



2. Basic scenario 

The basic scenario has been outlined in Fig. 1 of Araudo et al. (2007). A massive 
young stellar object (YSO), or a group of them, is embedded into a molecular 
cloud. The protostar heats the cloud in such a way that it can be detected as a 
strong IR source, with luminosities in the range LiR ~ 10^-^ Lq ~ 1038-39 

erg 

s~^ whereas the optical counterpart is obscured by the cloud. Masses and sizes 
of these clouds are of the order of ~ 10^ Mq and few pc, respectively (e.g. Garay 
&; Lizano 1999). The thermal jets can propagate up to distances of 10^^^^^ cm 
from the central source, and the non-thermal radio lobes are seen at distances 
of ~ pc from the protostar (Martf et al. 1993, Garay et al. 2003) with sizes 
of the order of 1 % this distance. The radio lobes are ionised by the terminal 
shocks and particles can be accelerated up to relativistic energies both in the 
reverse shock, as well as in the forward shock in the ionised gas. In what follows 
we will discuss the characteristics of the non-thermal particle population of the 
lobes arose from forward and reverse shocks. 

A crucial parameter for the acceleration of particles up to high energies is the 
magnetic field B. Equipartition calculations with the non-thermal population 
of particles (Araudo et al. 2007) and Zeeman measurements (e.g. Crutcher 
1999) suggest B ^ 1 mG. Magnetic fields are thought to play an important role 
supporting the cloud before the gravitational collapse, allowing high densities in 
the cores to be achieved (e.g. McKee &; Ostriker 2007). The medium density 
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strongly determines the jet evolution, and is a key factor regarding the emission 
at high energies. 

Here, we consider the effect of both a forward and a reverse shock, the 
escape of the particles from the acceleration region, and the illumination of the 
whole molecular cloud by protons that diffuse from the lobes, among others 
refinements with respect to our previous work. The cloud acts as reservoir of 
locally accelerated cosmic rays. In very massive clouds, where several protostars 
are forming, this effect can be very significative. 



3. Particle acceleration and losses 



Diffusive shock acceleration (DSA, e.g. Drury 1983) can operate in the turbulent 
region at the termination of the thermal jets. If the jet head is still moving in 
the free propagation phase, the bow shock in the cloud can be as fast as the jet, 
but then the shock luminosity must be smaller than the one carried by the jet 
(Lj). In any case, regardless the energy budget, the fast velocities allow particle 
acceleration up to very high energies. On the other hand, when the swept cloud 
material is massive enough, the jet slows down and the bow shock produced by 
the interaction can be as slow as ~ 10 km s~^. The shocked material of the cloud 
will be radiative, forming a very thin shell pushed by the jet. This shock is not 
suitable to accelerate particles, but its downstream material would be a good 
target for pp interactions or relativistic Bremsstrahlung in case that electrons 
and protons accelerated elsewhere could reach it. In the slow bow-shock case, 
the region suitable for particle acceleration is the reverse shock formed in the 
jet, which was not strong when the bow shock was moving as fast as the jet itself 
(being an effective accelerator). When the bow shock has slowed significantly 
down, the reverse shock has a velocity Vg ~ Ujet- The shock luminosity in the 
slow bow-shock case is similar to the jet one. 

In any case (bow shock or reverse shock acceleration), the particle acceler- 
ation rate for a parallel non-relativistic strong shock in the test particle approx- 
imation is 



dE, 



e,p 



dt 



rjqBc, 



gam 



20/s, 



(1) 



where e and p stand for electrons and protons, respectively (e.g. Protheroe 
1999). The parameter /sc is the ratio of the mean free path of the particles to 
their gyro-radius rg. Close to the Bohm limit, where the diffusion coefficient 
is Db = crg/3, /sc ^ 1. For outflows with bulk velocities v-.t^t ~ 10^ km s'^ 



'jet 

(e.g. Marti et al. 1995) we get a relatively low efficiency of r/ 10~' 
nonetheless, enough to achieve relativistic energies for the particles. 



which is, 



3.1. Hadronic losses 

The relevant energy loss for protons is due to pp collisions with the nuclei of the 
ambient gas of density n: 

= ncapp{Ep)KppEp, (2) 

pp 



dEj 
"dt 
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where app{Ep) = 34.3 + 1.88L + 0.251'^ mb, with L = ln{Ep/l TeV), is the 
interaction cross section (Kelner at al. 2006) and Kpp ~ 0.5 the inelasticity. 
The emission from pp interactions can be significant from high- to very high- 
energy 7 rays. 

In addition to the radiative losses, protons can escape from the emitting re- 
gion of size /lobe by diffusion or convection downstream the shock, in accordance 
with the timescales of these processes: tdiflf ~ ^lobc/^^B and tconv ~ 'ihohe/vs- 
The escape time tesc is defined as the minimum between tdis and tconv In the 
presence of a turbulent ~ 1 mG or higher, tconv will usually be shorter than 
tdiflf and tpp = Ep{dEp/dt)~p . Nevertheless, in the cases studied here the most 
energetic protons have their energy constrained by diffusive escape upstream. 

3.2. Leptonic losses 

The main losses for electrons are due to synchrotron radiation, IC up-scattering 
of IR photons, and relativistic Bremsstrahlung. The standard formulae are 
presented, for instance, in Blumenthal & Gould (1970). Synchrotron and IC 
losses in the Thomson regime satisfy: 



dE,, 



dt 



fdEe 



synch y 



V dt 



/ = ^>l, (3) 

IC/ ""ph 



where ub = -B^/Svr is the magnetic energy density (~ 10 ^ erg cm ^ for 
B = 1 mG), and Uph is the photon energy density inside the cloud, typically 
Uph ~ 10~^ erg cm~^, which is dominated by IR photons. This means that 
synchrotron losses will dominate over IC cooling. The highest energy electrons 
will cool via synchrotron and they will radiate from radio to soft X-rays through 
this mechanism. On the other hand, lower energy electrons are dominated by 
relativistic Bremsstrahlung and escape losses from the emitter. The relativistic 
Bremsstrahlung component contributes to the spectrum from hard X-rays up to 
high energies. 

For very high densities, ionization/ Coulomb collision losses may dominate 
radio emitting electrons, hardening their spectra. An observed steep radio spec- 
trum implies that this cooling channel is minor in comparison to synchrotron 
radiation. 

3.3. Maximum energies and particle spectra 

In the Bohm limit we expect for electrons that the maximum energy will be 
determined by the radiative losses. Equating energy gains and losses we get, for 
the parameters of the basic scenario sketched in Sect. [23 that E^^^ ~ 1 TeV. 
For protons, on the contrary, cooling times are longer and hence they can diffuse 
away from the shock upstream into the cloud. The residence time in the shock 
region, of size ~ hohe, determines the maximum energy: i?™^^ ~ 1 — 100 TeV. 

The particle distributions Ne^p{t, E^^p) are given by the transport equation 
(e.g. Ginzburg & Syrovatskii 1964): 

dNe,p {t,Ee,p) d[he,p{Ee,p)Ne,p {t,Ee,p)] Ne,p {t, Ee,p) ^ p ^ ... 
flT + + J = Qe,p{i^ ^e,p), (4) 
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where Qe,p{'t, E(, p) cx Ee,p'''^ is the function for the particle injection (a power 
law according to the DSA mechanism which is supposed to operate in the radio 
lobes). Injection takes place during the lifetime of the source (tiife), estimated 
typically in ~ 100 - 1000 yr « 3 x (10^ - 10^°) s for massive O protostars 
(see Garay et al. 2003). In the transport equation t is time, and be,p{Ee^p) is the 
summation of all energy losses. 

Eq. ([3]) can be solved in the steady state setting dN^^p (t, E^^p) /dt = 0. In 
the cases when synchrotron dominates over escape or relativistic Bremsstrahlung, 
the electron distribution will be Ne{Ee) oc QeE~^ at high energies. Relativistic 
Bremsstrahlung losses and escape, on the contrary, lead to Ne{Ee) oc Q^{Ef,). 
The injected proton spectrum is not modified since the relevant losses, due to 
inelastic pp collisions, are almost linear with the energy in the range considered 
(from ~ 100 MeV to TeV) Np{Ep) oc Qp. 

The normalization of the electron spectrum is given by the observed radio 
luminosities, between 10^*^ and 10^*^ erg s~^ for sources like HH 80-81 (Marti et al. 
1993) and IRAS 16547-4247 (Rodriguez et al. 2005). The proton to electron 
energy density ratio, a = Up/ue, cannot be established from first principles 
in the acceleration process and it should be considered as a phenomenological 
parameter. All this fixes the luminosity injected in non-thermal particles (Lnt)- 
In the next section we will consider two cases: one with only electrons (a = 0) 
and another dominated by protons, with a ratio a = 100, as observed in galactic 
cosmic rays. 



4. Non-thermal spectral energy distributions 

Once the particle energy distribution N(,^p{Ee^p) has been obtained through Eq. 
dl]), the spectral energy distribution (SED) can be calculated as: 

E-jLe^ = E'^ j qj{E^, Ee,p)Ne^p{Ee^p)dEe^p, (5) 

where is the particle emissivity of the particular radiation process and Le^ 
is the specific luminosity. 

4.1. IRAS 16547-4247 

The powerful IR source IRAS 16547—4247 is composed by a massive protostar 
embedded in a very dense and luminous molecular cloud (Garay et al. 2003). In 
addition, two radio lobes are at the end of the thermal jets of the central YSO, 
being the south one clearly non-thermal. 

In Fig. [1] (Left) we show an example of the SED for the non-thermal lobe 
of the IRAS source, where only primary electrons have been considered. In this 
specific case, the reverse shock is accelerating particles since the cloud material 
is very dense and the bow shock has been already slowed down. Densities have 
been taken as high as possible accounting for the ionization cooling limit. The 
parameters used in the modeling are presented in Tabled) This example shows 
that, despite the IC contribution is negligible, 7-rays might still be of leptonic 
origin if densities are particularly high. The spectral shape is determined by 
the escape of electrons from the emitter (this was not considered in Araudo 
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et al. 2007). The SED has two clear peaks, one due to synchrotron radiation 
and the second one due to relativistic Bremsstrahlung. The latter emission 
extends into the 7-ray domain with integrated luminosities of ~ 5 x 10^^ and 
3 X 10^^ erg s"^ at ~ 0.1 - 100 and > 100 GeV, respectively. The emission 
for the case with a = 100 has been also estimated but not shown in the present 
paper. Synchrotron emission is now dominated by secondary pairs, and the pp 
in the lobe reaches integrated luminosities of ~ 2 X 10^3 and 2 x lO^^ gj.g s"^ at 
~ 0.1 - 100 and > 100 GeV, respectively. 

4.2. HH 80-81 system 

The Herbig-Haro objects HH 80 and HH 81 are located at the terminal point 
of a massive YSO jet. The triple system composed by the central YSO plus 
HH 80-81 and HH 80 North is at the edge of a dark cloud (Marti et al. 1993). 

In Fig. [1] (Right), we show a case dominated by protons, with a = 100, 
for the source HH 80. Given the lower densities of the environment, the bow 
shock could still move fast and therefore be an efficient particle accelerator. In 
this source the radio emission is still produced by primary electrons, with some 
contribution from secondary pairs injected through the decay of charged pions. 
At high energies the IC and the relativistic Bremsstrahlung are negligible and 
the emission is dominated by 7 rays originated in neutral pion decays from pp 
interactions. The 7-ray cutoff appears at energies E^^^ ~ The total 

7-ray integrated luminosities from the lobe can reach ~ lO'^'^ erg s~^. However, 
since protons can diffuse out the acceleration region with significant energy due 
to their long cooling times, the interaction with the cold material of the cloud will 
add a diffuse component (Aharonian & Atoyan 1996). This emission could look 
marginally extended for current Cherenkov telescopes and Fermi. The diffuse 
radiation from the massive molecular cloud can increase the 7-ray integrated 
luminosity to the level of ~ 5 x 10^^ and 6 x 10^^ erg s~^ at E'-y ~ 0.1 — 100 
and > 100 GeV, respectively, rendering the source detectable. The leptonic case 
(a=0) for the object HH 80 was also computed but the luminosities at high 
energies are well below lO^'^ erg and the spectra are not shown here. 

Table 1. Physical parameters considered in this work. 





IRAS 16547-4247 


HH 80 


Cloud 


^lobo [cm] 


1.1 X 10^*^ 


5 X 10^^ 




Lj [erg s"i] 




1.5 X 10^^ 




tdyn" [S] 


4 X 10^ 


3 X 10^ 


lOii 


a 





100 




n [cm~3] 


10^ 


10^ 


2.5 X 10^ 


B [mG] 


1 


0.3 




particle 


e 


P 




0.02 


0.1 




r 


2.2 


1.8 





" Dominant timescale: it is approximately tiife for HH 80 and the cloud, and t, 
for the IRAS source. 
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Figure 1. Left: SED of the south lobe of the source IRAS 16547-4247 
for a pure leptonic case (a=0). The IC contribution is negUgible and not 
shown here. Radio and X-ray observational points are from Rodriguez et 
al. (2005) and from Araudo et al. (2007), respectively. Right: SED of the 
source HH 80 for a case dominated by protons (a = 100). The synchrotron 
radiation is dominated by primary electrons. Radio and X-ray observational 
points are from Marti et al. (1993) and Pravdo, Tsuboi & Maeda (2004), 
respectively. Both the emission from the lobe and that which results from 
the proton diffusion in the cloud are shown. In both figures we indicate the 
sensitivity curve of Fermi for 1 yr of observations and the distance of the 
sources, - 3 and 1.7 kpc for IRAS 16547-4247 and HH 80-81, respectively. 



5. Summary and discussion 

In the present contribution, we study the production of 7-ray emission at the 
terminal points of jets that emanate from massive protostars. In particular, 
we consider two sources: IRAS 16547-4247 and HH 80-81. In the former case, 
ralativistic particles are accelerated in the shocked jet, wereas in the latter, 
in the bow shock. With the model presented here, 7-ray luminosities up to 
~ lO'^'^ erg s~^ are obtained. 

As pointed out by Romero (2008), the detection of massive protostars at 7- 
ray energies would open a new window to star formation studies. The detection 
of the cutoff in the SED will give important insights on the acceleration efficiency 
in the terminal shocks of the outflows. This, in turn, can be used to constrain 
the shock velocity and the diffusion coefficient. The accumulation of cosmic rays 
accelerated in the radio lobes into the molecular cloud can produce extended 7- 
ray sources. 

The combined effect of several protostars deeply embedded in giant clouds 
might be responsible for GeV-TeV sources found in star forming regions by 
EGRET, Fermi, AGILE and Cherenkov telescopes. A clustering of gamma-ray 
sources should be present in regions with large molecular clouds and star forma- 
tion, as already inferred from EGRET data (e.g. Romero et al. 1999). Although 
we do not expect that massive protostars should be among the bright sources 
detected by Fermi (Abdo et al. 2009), our predictions show that they should 
show up in further analysis of weaker sources after few years of observations. 
The emission levels above 100 GeV, of about 0.01 Crab, could be detectable by 
current and future Cherenkov telescopes for ~ 50 hours observation time. 

Cosmic-ray re-acceleration inside the clouds due to magnetic turbulence 
(e.g. Dogiel et al. 2004) could result in stronger sources. Neither UV nor hard 
X-ray counterparts related to thermal Bremsstrahlung produced in the shock 
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downstream regions are expected to be observed from these sources because 
of the large absorption and/or low emission levels. In any case, massive clouds 
with high IR luminosities and maser emission (tracers of massive star formation) 
deserve detailed study with Fermi and other 7-ray telescopes. 
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